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Abstract 

We propose a new type of magnetically disordered ground state for a frus- 
trated quantum antiferromagnet. This disordered state is an array of spin 
singlets spontaneously formed on four spin plaquettes. Both perturbation re- 
sults and bond-operator calculations show that this phase has lower energy 
than the columnar dimer state. Analysis of available numerical data on finite 
clusters also supports the conclusion that this state is realized at intermediate 
frustrations. 

PACS numbers: 75.10.Jm, 75.10.-b, 75.30.Kz, 75.50.Ee 



Typeset using REVT^ 



Frustrated two-dimensional quantum magnets are a fascinating topic of numerous studies 
over the past decade. The interest has primarily focused on melting of long-range magnetic 
order and appearance of disordered phases at T = as a result of enhanced quantum 
fluctuations. The most popular spin system having such type of behavior is a frustrated 
square lattice Heisenberg model with the nearest-neighbor antiferromagnetic exchange Ji 
and the second-neighbor coupling J2 

H = Ji5:Si-s, + JsEs.-s, . (1) 
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Analysis within the linear spin-wave theory reveals that for weak frustration the model has 
a Neel ground state with ordering wave- vector (7r,7r), while for strong frustration spins are 
ordered at wave- vector (7r,0) or (0, tt) (stripe or coUinear state). For any finite S there is 
also a region around the classical critical point J2/ Ji = 0.5 where sublattice magnetization 
vanishes This result raises a question about the ground state for intermediate frustrations 
which may be nonmagnetic and represent a kind of two-dimensional spin liquid. 

Subsequent theoretical works have considered the problem of the intermediate phase of 
the frustrated antiferromagnet Eq. (|I|) mainly from three different points of view. The first 
group includes exact diagonalization studies on finite clusters, which clearly show existence 
of a disordered spin state in the region 0.4 < J2/J1 < 0.65, though final conclusion about 
its nature has not been reached The second group of works address the problem 

by calculating higher-order 1/S corrections either in framework of the modified spin- wave 
theory or using the Schwinger boson mean- field calculations Both approaches predict 

enhanced stabilities for the Neel and the coUinear states resulting in a finite overlap of the 
magnetically ordered phases. Second-order corrections to the mean-field solution, on the 
other hand, suggest a small window 0.52 < J2I J\ < 0.57 between the two ordered states for 

The third group of works explore one particular possibility for disordered ground state — 
a valence bond crystal with broken translational symmetry and a finite gap in excitation 
spectrum. Studying nearest-neighbor SU(A^) antiferromagnet in the large limit. Read 
and Sachdev have argued in favor of a state with spontaneous columnar dimerization [0. 
Subsequent calculations for the spin-i Hamiltonian (|I]) using series expansion [§] and boson 
techniques supported stability of this phase around J2/J1 = 0.5. Consistent with these 
treatments cluster results for dimerized susceptibility yield its noticeable increase in the 
intermediate region @J^. 

In this Letter we propose a new type of valence bond crystal ground state for a frustrated 
quantum antiferromagnet. This is a so called plaquette resonating- valence-bond (RVB) state 
shown in Fig. 1(a). We calculate its energy by two methods and find that it is lower than 
for the columnar dimerization [Fig. 1(b)] in the same approximation. Then, by considering 
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relevant order parameters, we show that the existing numerical data on the 6x6 cluster 
are, as a matter of fact, in favor of this state. 

Two types of magnetically disordered singlet ground states have been widely discussed 
for a spin- 1 square lattice antiferromagnet. The first, a featureless RVB spin-liquid, describes 
a linear superposition of spin pair singlets with long- or short-range correlations [0 . The 
second proposal is a spin-Pierls order of valence bonds which are frozen and break lattice 
symmetries 0. One can also consider an intermediate situation: valence bonds resonate in 
finite spin blocks resulting again in breaking of translational symmetry. In the simplest case 
such spin blocks contain four sites. Similarly to discussion of dimerized phases, the first step 
is to investigate a restricted Hamiltonian of four spins on a square plaquette: 

Hplq= Jl(Sl + S3)-(S2 + S4) + J2(Sr 83 + 82- 84). (2) 

This Hamiltonian is easily diagonalized. The ground level for J2 < Ji is a singlet character- 
ized by the quantum numbers S'13 = 5*24 = 1 with the energy Eg = —2Ji + ^J2- Its wave 
function is \s) = -^({1, 2}{4, 3} + {1,4}{2,3}), where curly brackets denote singlet bond 
of a spin pair. This state can be considered as a RVB-like state in a four spin subsystem 
with energy lower than for a frozen dimer configuration, e.g., {1,2}{4, 3}. Excited levels 
are three triplets \ta), \Pa), IQo) = x,y,z), Et = —Ji + \ J2, Ep = Eg = —\J2, a quintet 
\du) (u = 1,...,5), Ed = Ji + |J2, and another singlet \s') = {1,2}{4,3} - {1,4}{2,3}, 
Eg' = — 1^2, which crosses with |s) and becomes the ground state for J2 > Ji (we will not 
be concerned with this case in our paper). 

Four spin singlets form a plaquette covering of a two-dimensional lattice. Previously, 
the plaquette-RVB state has been predicted for the disordered two-dimensional magnet 
CaV409, which is described by a spin-| Heisenberg model on a 1/5-depleted square lattice 
|[TT| . In that case singlets correspond to spin blocks chosen by the lattice geometry. For 
the translationally invariant Hamiltonian ([l|) periodic array of plaquette singlets appears 
spontaneously. The ground state has in this case four-fold degeneracy determined by broken 
translations along two sides and diagonal of an elementary square. Using simple rules for 



products of dimer coverings [T^ one can show that the overlap between the two states in 
Fig. 1 formed by pure singlets decreases to zero as (a/3/2)^/^ with increasing number of 
spins A^. Consequently, the plaquette-RVB state cannot be represented as a superposition 
of two column states rotated by 90° with respect to each other. Later, we consider spin 
order parameters, which have different values for the two disordered phases. 

Since the Hamiltonian (|l|) includes interaction of spins from different blocks, the ground 
state does not coincide with a simple product of block singlets and nonzero expectation val- 
ues (Sj ■ Sj) appear between all nearest neighbor pairs. To compare energies of the columnar 
dimer and the plaquette-RVB phases we use a boson technique suited for perturbative ana- 
lytical expansion around the local spin singlets. It was used previously to study dimerized 
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phases [||. We first generalize this method by deriving simultaneously boson representations 
for dimer and plaquette spins and then calculate energies of the two states in the mean-field 
approximation. 

Let us consider an arbitrary spin block with singlet ground state and excited levels 
denoted by A spin in a block is expressed in terms of the basis states by 

where Z^^ is the projection operator n is a local spin index inside the block, and 

the global block index in the lattice i is omitted for simplicity. We derive first the matrix 
elements in a subspace of the ground state singlet and the lowest triplet states. These 
states form a complete set for dimers, whereas for plaquettes this is only a part of the local 
basis. Later, we partially take into account higher energy states of plaquettes. Rotational 
invariance in spin space and time-reversal symmetry give 

(tQ,|S^|t^) = %e ^'^ A^^ , 

where e°'^"' is the totally antisymmetric tensor and A"^, A"^ are real constants. Us- 
ing the explicit forms of singlet and triplet wave functions one straightforwardly gets 
A^^ = (-l)"/2, A'l^ = 1/2 for dimers, and = (-1)"/V6, = 1/4 for plaquettes. 
Let us now define the vacuum |0) and four boson operators which yield the four physi- 
cal states by \s) = s~^\0), \ta) = t^\0). The projection operators are, then, expressed as 
^sta — s+t^^ ^tatfj _ t+tp^ and so on. Block spins represented via these boson operators are 

= ^{sn^ + tts)- ie"^^t+t^, for dimers, 

= i^{s+ta + tts)- {e^^^tp^, for plaquettes 

Spin commutation relations are satisfied in the chosen subspace as long as the bosonic 
representation preserves the algebra of the projection operators: Z^^ Z^'^' = S^iyZ^^' . This 
requirement restricts the boson occupation numbers: 

s^s + Y.tlt^ = l. (4) 

a 

With this constraint the Hamiltonian of a single block, a dimer or a plaquette, takes a form, 
Tis = EsS~^s + Ettp^. The simplest way to deal with Eq. @ is to impose it as 'a constraint 
in average' enforced by a chemical potential. 

Further calculations follow closely the work by Sachdev and Bhatt |p. Using Eq. we 
rewrite the Hamiltonian (|l|) in terms of bond operators for each of the two states and assume 
site independent chemical potential /i and condensate of singlets (sj) = s. Interaction terms 
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in the boson Hamiltonian are classified by the number of triplet operators t^k- It can be 
shown that the terms with three and four triplet operators affect the results only slightly 
1^. Therefore, we omit them and after diagonalization of the remaining quadratic form 
determine the parameters fi and s through the saddle-point equations 

dEg,,Jdfi = , <9Eg.s./(9s = . (5) 

The system is a magnetically disordered phase as long as the gap in the excitation spectrum 
is positive. Vanishing of the gap leads to a condensation of triplet excitations: (tjQ,) 7^ 0. 
Note, that this does not mean transition, e.g., to the usual Neel state, since spontaneous 
dimerization of bonds will be preserved in such an ordered phase. 

For the columnar dimer state we recover the results of 0. This state is stable in the 
region 0.19 < J2/ Ji < 0.66. For the plaquette-RVB phase analogous treatment would under- 
estimate the energy because of the neglection of the higher-lying states of plaquette. Among 
the other excited levels the most significant energy corrections are given by higher triplets 
which have nonzero matrix elements with the ground state for spin operators. Extending 
calculations preceding Eq. to these excited triplets we find the following representation 
of spins in each plaquette 

(6) 

-^2,4 = -75(^^^« + its) ± + q^s) , 

with the constraint (|D being changed to s~^s + J2a(tt'l^a + ptiPa + itla) = 1- According to 
the quadratic approximation we have omitted in Eq. @ products of two triplet operators. 
The boson Hamiltonian obtained with this substitution can be further simplified by keeping 
interaction terms with only one of the higher triplets pa and Qa- Physically, it means that 
we take account of higher-lying triplets only via scattering of lowest excitations on them. 
As a result, one mode in the upper band remains constant, while the other acquires a finite 
dispersion. 

Solution of the saddle-point equations (j|) gives us parameters ^ and s, from which 
we calculate energy and spin gap of the plaquette-RVB phase. This disordered state is 
locally stable for J2 > O.O8J1. For large frustrations stability boundary of this phase lies 
at J2 > O.8J1, well inside the region where the coUinear state is expected to appear. This 
boundary cannot be determined correctly in our approximation as it depends crucially on 
the behavior of the neglected higher singlet. Therefore, we restrict ourselves to the region 
J2 < 0.7Ji. The gap of the plaquette state at J2 = 0.5Ji is Apiq = 0.85Ji, while for the 
columnar dimer state Admr = 0.74Ji. The results for the energies of these two phases in 
units of Ji are presented in Fig. 2 by solid lines. In addition to significantly wider region of 
stability, the plaquette-RVB phase has lower energy than the dimerized state for J2 < 0.58 Ji. 
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At J2 = O.5J1 the difference between the two energies, Ep\q = —0.466 and i?dmr = —0.456, 
is about 2%. (Without higher triplets Ep\q = —0.458.) Thus, the bond-operator formahsm 
predicts the plaquette-RVB state instead of the columnar dimer state as the intermediate 
phase of the frustrated antiferromagnet (|1|). 

Another way to estimate analytically the energies of the two valence bond states is 
the second order perturbation expansion starting from the singlets either on dimers or on 
plaquettes. This approximation corresponds to the first two nonzero terms in the series 
expansion method p|. The results are shown in Fig. 2 by dashed lines. At J2 = 0.5Ji 
the difference between E^i^ = —0.63 and -Edmr = —0.492 is much bigger than in the bond- 
operator scheme. Though significant corrections are expected in higher orders of series 
expansions, these calculations agree with our proposal of the plaquette-RVB state for the 
ground state of the Hamiltonian (|I|). They also show that the crossing of the two valence 
bond levels found by the boson technique may be shifted to stronger frustrations. 

Both results for the energies of the two disordered phases have been obtained by ap- 
proximate analytical methods and, hence, may be questioned. We, therefore, now discuss 
how the columnar and plaquette phases can be distinguished in exact numerical diagonal- 
ization studies. The key step is to construct appropriate order parameters, which measure 
nonequivalence of spin bonds in a disordered state. A set of such spin operators have been 
proposed by Sachdev |]12|. Appearance of the column phase has been tested 0,0] by using 
the parameter 

This operator is site independent and takes values proportional to 1, i, —1, and —i for 
the four degenerate column states. On the other hand, for the four plaquette-RVB states 
"^i equals e*'^/'^, e*^'^/^, — e*'^/^, and — e*^'^/'^ up to a constant factor. Since in numerical 
studies on finite systems one measures only the corresponding correlation function Xtp = 
1/N'^{\ I^i^l'iP), it is impossible to distinguish the two states by observing an anomaly in 
XipiJ^)- To overcome this difficulty another order parameter should be considered 

$i = Si ■ (Si + y + Si_y - Si + x - Si_5c) . (8) 

For the four column states $i has values 1,-1,1, and —1, whereas it becomes zero for each of 
the plaquette-RVB states. Consequently, the two crystalline bond states can be distinguished 
by measuring two quantaties Xtp and x</> simultaneously. The order parameter (||) coincides, 
in fact, with the magnetic structure factor at the wave vector (vr, 0) and has been studied 
to check the difference between the two ordered magnetic phases: $i = in the Neel state, 
$i 7^ in the collinear state. Thus, numerical results for both and X(t> currently 
available. We refer to the most reliable data obtained on the 6x6 cluster ||^. Fig. 3 of this 
work plots two correlation functions Xi(>{J2) and Xii>{J2)- It is clearly seen that at the peak 
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point of Xi>{J2) the other function differs only shghtly from its value at J2 = 0, whereas 
X<j>{J2) starts to increase simultaneously with decreasing of Xi(>- This observation is a direct 
convincing evidence of the plaquette-RVB ground state for the model (|I|) at intermediate 
frustrations. 

We now turn to physical consequences of the translational symmetry breaking in a mag- 
netically disordered ground state of spin system. The first property we address is a nature 
of the phase transitions at T = 0. If transformation from such a disordered state to a 
spin ordered state, e.g., to the Neel state, is continuous, the usual group-subgroup relations 
should be satisfied for symmetry groups of the two phases. However, this is not the case for 
the Neel and the plaquette-RVB (or columnar dimer) states, since neither the former has 
a higher symmetry than the latter (because of the absence of rotational invariance in spin 
space), nor the latter is more symmetric than the former (due to translational symmetry 
breaking). Analogous arguments are applied to the transition into the coUinear state. Thus, 
the magnetically disordered ground state region of the model (|l|) is bounded by two points 
of first order transitions. This conclusion explains, in particular, diverse estimates of its 
width @-§|,PJ^ found by studying stabilities of the different phases. 

Finally, we comment on finite temperature behavior of valence bond crystals. As they 
break only the discrete lattice symmetries, there should be transition from the symmetric 
paramagnetic phase with decreasing temperature. Analogous conclusion for the collinear 
state has been reached in |T^. The difference between these two cases lies in symmetry 
properties of the corresponding order parameters. At T 7^ 0, when sublattice magnetization 
of the collinear phase vanishes, the order parameter is a soft Ising type quantity $i having 
values +1 and —1 on each lattice site. In this case only the rotational lattice symmetry is 
broken in the low-temperature phase. Singlet formation in the bond crystalline phases is 
accompanied by translational symmetry breaking and corresponds to k 7^ 0. Bond strength 
modulations appear at wave vectors (vr, 0) or (0, vr) for the columnar dimerization and the 
spin bond order parameter has two components: Si • (Si_x — Si+x) and Si • (Si_y — Si+y). For 
each of the four column phases only one component is nonzero. The plaquette-RVB state 
also corresponds to these wave-vectors having both components nonzero at the same time. 
One can straightforwardly construct a two-component Landau free-energy functional for this 
irreducible representation and show that the choice between (1,0) and (1, 1) symmetries is 
determined by minimization of the 'anisotropic' fourth order term. Our calculations for the 
energies of these two states at T = suggest that the minimum occurs for (±1, ±1) states. 

We thank O. A. Starykh and M. Troyer for helpful discussions. This work has been fi- 
nancially supported by a Grant-in-Aid from the Ministry of Education, Science and Cultutre 
of Japan. 
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FIGURES 

FIG. 1. (a) plaquette-RVB and (b) columnar dimer ground states; bold lines denote stronger 
spin bonds. 

FIG. 2. Energies (in units of Ji) of (a) columnar dimerized and (b) plaquette-RVB phases 
calculated by bond-operator technique; (c) and (d) for the same phases by the second-order per- 
turbation theory. 



9 



(a) 



(b) 



